A Field Theory for Finite Dimensional 
Site Disordered Spin Systems 
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We present a new field theoretic approach for finite dimen- 
sional site disordered spin systems by introducing the notion 
of grand canonical disorder, where the number of spins in the 
system is random but quenched. We perform the simplest 
non-trivial analysis of this field theory by using the variational 
replica formalism. We explicitly discuss a three dimensional 
RKKY-like system where we find a spin glass phase with con- 
tinuous replica symmetry breaking. 

PACS numbers: 75.10N, 75.30F, 75.50L 



Most advances in the field of disordered spin systems 
have been based on models in which the bonds take ran- 
dom values. However in most experimental systems the 
positions of the spins are random but the interactions oc- 
cur through deterministic potentials. Analytic studies of 
site-disordered spin systems, such as RKKY spin glasses 
and dilute ferromagnets, have been hampered by the lack 
of a suitable field theoretic model (however a lattice based 
formulation has been proposed M). By considering a 
situation in which the number of spins in the system 
is random but quenched we are able to write a replica 
field theory for site-disordered systems. This field the- 
ory seems to be simpler than many of those coming from 
bond-disordered and diluted lattice models and should be 
accessible to many standard analytical techniques. The 
mean field theory of this model, for reasons that will be- 
come clear, cannot provide any information about spin 
glass order. In the second part of this letter we con- 
sider the simplest generalisation of mean field theory, the 
Gaussian variational (GV) method, which does provide 
this information. Use of the GV method is widespread, 
and it is a useful warning that for certain interaction 
types in our model it gives unphysical predictions. 

The role of replica symmetry breaking (RSB) in disor- 
dered spin systems is of great interest. Although RSB in 
the mean field theory for spin glasses is now well under- 
stood ||] and related to the proliferation of pure states of 
the system, in finite dimensions the picture is less clear. 
Alternative qualitative approaches based on droplets H 
view the spin glass phase as a disguised ferromagnetic 
phase with only two underlying fundamental states. We 
will explicitly consider a 3-dimensional site disordered 
spin glass using an RKKY-like interaction in our model 
and find continuous RSB in the GV approximation. 



Although in this letter we concentrate our attention on 
spin glass physics with oscillating sign interactions, the 
model can also describe dilute ferromagnetic or antiferro- 
magnetic systems. Indeed, even for the RKKY example 
we find ferromagnetic order at very low temperatures. 
The application of the methods described here to these 
single sign interactions is an interesting subject, but we 
defer it to a longer article Q , simply mentioning some of 
the issues that arise at the end of this letter. 

Firstly consider a model where the number of spins N 
is fixed: N spins Si are placed randomly at positions 
uniformly throughout a volume V. This type of disorder 
we refer to as canonical disorder, as the number of parti- 
cles is the same for each realization of the disorder. The 
spins interact via a pairwise potential J depending only 
on the distance between the spins. The Hamiltonian is 
then given by 



Tj)SiSj 



(1) 



Assuming that J is positive definite, making a Hubbard- 
Stratonovich transformation expresses the partition func- 



tion as 



(2) 
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Employing replicas, we average out the site-disorder by 
integrating over the positions rj using the flat measure: 

Z n N = JV(t> a [detJl3]-2 x 

expf-^ / [j^M^J-Hr-r'W^drdr' 



Vlogi J f[ 2 cosher) dr 



(3) 



A field theoretic analysis of the above theory is compli- 
cated by the presence of the log term in the action. We 
overcome this difficulty by making a physically desirable 
modification to the definition of the disorder. In general 
one might expect the system to have been taken from a 
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much larger system with a mean concentration of spins 
per unit volume, p. A suitably large subsystem of volume 
V will thus contain a number of spins N which is random 
and Poisson distributed: p(N) = exp(— pV)(pV) N /N\. 
This distribution must be used to weight the averaged 
free energy so we are led to define H n = J2nP(N)Z n . 
By analogy with the statistical mechanics of pure sys- 
tems, we shall call this type of disorder "grand canonical 
disorder" . 

The resulting theory is simpler than (j^) and is defined 

by 

S" = exp(-pV) J V(t> a {detJ(3}-% x 

ex P (-^ J jY,Mr)J- l {r-r')4> a {r')drdr' 



][[ 2 cosher) dr). (4) 



Expanding the cosh one sees that the leading term cor- 
responds to the random temperature or random mass, 
familiar from bond disordered approaches, and that de- 
pending on the choice of interaction one might expect 
similar renormalisation group results QJ^] . 

In order to relate this theory to measurable quantities 
we return to the original formulation of the model in 
equation and identify physical operators. The spin 
density operator, M a (r) — S(r—ri)Sf is closely related 
to the field 4> a appearing in the theory. The equations of 
motion following from the replicated version of (^|) show 
that the physical magnetisation density is given by, 



1 r i N 

M =- jdr{M a {r)) = -{Y,St) = 



f3J(0y 



(5) 



and the correlator (M a (r)Mb(r')) , is in terms of G a b(r - 
r') = {<f> a {r)<j>b(r'))c, obtained as 



(M a (k)M b (-k)) c = 



gab 



G ab (k) 

(3 2 J 2 (k) ~ f3J(k) 



(6) 



M a (r) is not, however, the operator sensitive to spin glass 
ordering, and it is natural to consider another operator 
Qab{ r ) = J2i 3{ r ~ r a)S^S b , related to the non-linear sus- 
ceptibility. This new operator is composite and does not 
manifestly appear in the field theory (Q), it is for this rea- 
son that we must go beyond mean field theory to obtain 
non-trivial results. Operators involving more spins can 
be introduced in the same way. 

In the remainder of the letter we analyse this field the- 
ory with the Gaussian variational method which can be 
regarded as a generalisation of mean field theory. This 
method, otherwise known as Hartree Fock, is a trunca- 
tion of the Schwinger Dyson equations and becomes ex- 
act in the limit of many spin components (such an m- 
component theory is treated in a separate publication 



M). In the context of disordered systems, this method 
has has success in calculating exponents for random man- 
ifolds 0, but one should bear in mind that important 
effects may occur at higher orders in 1/m. In fact, for 
certain choices of the potential in the field theory consid- 
ered here, one can rigorously demonstrate a failure of the 
method 0. We return to a discussion of the reliability 
of the approximation at the end of the letter. 

We allow the possibility of ferromagnetic order and 
make the ansatze that (</> a (r)) = (j> a and {4> a {r)4>b{r')) c = 
G a b{r — r') (by translational invariance). The variational 
free energy is given, up to constant terms, by 



n(3F var = -^TrlogGU + -?-Tr G^J" 1 



(7) 



where the above traces are both functional and on replica 
indices and where f2 is defined by 



^exp ($^«£« + i$ZG o6 (0)S a S , 6] 

S a \ a ab ) 



(8) 



As usual we do not expect breaking of replica symme- 
try on single-index objects and hence set <j) a = 4>. The 
variational equations are 



and 



4> a = $ = i3 P J(o)n a 



Gab( k ) = -3^abJ- l {k)- P n ab . 



(9) 



(10) 



where fi a and £l a b are traces of the type (g) containing 
respectively S a and S a Sb- 

Within this approximation, by introducing a source 
for the operator qabir) and using the FDT theorem, we 
obtain an equation for the correlation function Q a bcd — 
(q a b(r)q c d(0))- 

^6 abed 

(k) = pflabcd + o ^abgh{k)Qghcd{k) 
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/d d p 
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where fl a bcd is another object of the type (^]) containing 
four S's. 

A replica symmetric (RS) ansatz for G leads to a 
regime specified by two order parameters: the magneti- 
sation M (jfy and the Edwards Anderson order parame- 
ter qEA- These parameters are determined by a pair of 
equations very similar to the mean field equations for the 
Sherrington Kirkpatrick (SK) model J| 
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q = 



^e-^-tanh 2 (/3J(0)M + £,^fgl) ) ( 1. 2 ) 



Where g\ (the off diagonal part of G a &(0)) is given by 
d d k J 2 {k) 



9i = pP q 



{2n)d (l - (1 - q) P PJ(k) 



(13) 



The simplest solution of these equations yields the 
high temperature, low density paramagnetic region with 
qEA — and M = 0. In this region, the two index cor- 
relator (^) is known and equation ( |TT| ) can be solved for 
the most interesting correlator, (q a b(k)q a b(— k)) 



pP 2 



1-p/ 



d d p 



J(k)J(p-k) 



(14) 



( 2 « d (l~(l-?)p/3J(fe))(l-(l-g)p/3J(p-fe)) 



This correlator is simply related to (S(k)S(-k)) 2 = 
12' a b(q ab (k)qab(— k)) and the divergence in the above for- 
mula signals the onset of a spin glass phase. The diver- 
gence occurs on a line in the temperature density plane 
specified by g± — q and coincides with the AT line as de- 
termined by stability considerations [[|||. Furthermore 
the phase boundary also coincides with the line on which 
the RS equations ( |l2| ) develop solutions with non-zero q. 
This situation also occurs in the SK model and suggests 
a continuous breaking of replica symmetry. 

We shall look for continuous replica symmetry bro- 
ken solutions and parameterise the off-diagonal part of 
the matrix G a b(k) by a continuous Parisi function g(k, u) 
where u G [0, 1], and a diagonal part denoted by gu(fc). 
For such a matrix, fl (J|) is very similar to the free energy 
in the SK model; it cannot be obtained in a closed form 
and a standard strategy is to work close to the transition 
line by expanding f2 up to a term of 0(g i ) which in the 
SK model is the first term leading to a breaking of replica 
symmetry. The expansion is JHJ 



n-i 1 



9d(0) 



(15) 



The remaining terms in the action are easily computed 
within the algebra of Parisi matrices JtJ . The variational 
equations one obtains are 



[g D (k)}- 1 = +pa D (k)= P r^—-l 



ftp 



(16) 



Defining 



D D (k) =a D + (a) 



D(u,k) =a D + (a) + [cr](u), 



(17) 



(in the notation of [Q), the equations can be inverted to 
find 

9d{k) - pDD ( k ) \ y -I- Jo ipr D (u,k) + D D (k) ) 

g(u,k) = p£) ^( fe ) (uD(ui) +/o IF D(l,k) + Do(k)) ■ 

Proceeding by differentiating the second equation of ( |l6| ) 
with respect to u one obtains a' = or 



1+5 -ug- g) = p 



d d k 1 
(2^7 D 2 



(19) 



Taking a second derivative in some region where equation 
(H) holds we find: 



/ 



a(u)u = — 
y ' 2 



1 + 2^ 



I 



d a k 



( r d d k 1 
\J (2ir) d TP 



(20) 



For 4 or more dimensions, in the limit in which the short 
distance cutoff is removed, this equation is simple and 
we find a scenario similar to that found in the SK model. 
In general the function a(u) depends on a{u) and one 
obtains a first order nonlinear differential equation for 
g(0,u). In all cases we have considered, the power se- 
ries solution near the origin starts with a linear term. 
For consistency with our perturbative analysis the region 
where g(0, u) is non constant must be close to the origin. 
More precisely there must be a break point u with small 
value, above which g(0,u) is constant and equal to au 
if the solution is to be continuous. The breaking pattern 
scenario is reminiscent of the random manifold problem 
with long range disorder and is qualitatively the same 
as found in the SK model near T c Q . 

It is useful to illustrate these results for a specific in- 
teraction and for the purposes of this letter we consider 
an RKKY-likc oscillatory potential in 3-dimensions. 



J(k) = p- 3 6(p-\k\). 



(21) 



The dimensional constant p merely sets the scale of the 
problem and can be set to 1. Using equation (|TJ) we 
obtain for (q a b(k)q a b(—k)} in the paramagnetic phase 



labab 



(d73j) (fe+4)(fc-2) 2 ' 



for k < 2; 
for k > 2. 



(22) 



The spin glass phase boundary is given by p = (1 + 
12tt 2 T- Vl + 247r 2 T)/(127r 2 ), and the exponent associ- 
ated with the transition is 77 = 1. Numerical inspection 
of the RS equations finds stable ferromagnetic solutions 
at low temperature because at very high densities the 
positive short range part of the potential can dominate. 
We illustrate the expected form of the phase diagram in 
figure 1. 



3 



The expansion just below the AT line gives rise to a 
differential equation as described above, the leading so- 
lution at small u being linear. The break point uq can 
be calculated in terms of the deviation from the AT line: 
i*o = 5/3^2 (i +p p)(i- p p) ■ Despite having the full struc- 
ture of the two index correlators g(k,u), leading to non 
trivial momentum dependence in related to the 

connected magnetic correlation function (Q) , the analysis 
only holds close to the spin glass transition and is un- 
able to address the ferromagnetic transition which takes 
place at much lower temperature. The four index cor- 
relation functions Q a bcd contain much of the physics of 
the spin glass phase: for example the 0-exponent || may 
be extracted from the long distance behaviour of such 
objects. Equation ( |il"| ) is however an equation carrying 
four replica indices and the solution in the case of con- 
tinuous replica symmetry breaking is technically rather 
formidable requiring extensions of the methods described 
in g§. 

As we have emphasised, we have used the Gaussian 
variational method because it is the simplest generalisa- 
tion of mean field theory that gives us access to spin glass 
physics. We now discuss the reliability of the approxi- 
mation. Certainly we should expect it to be exact for 
m-component Heisenberg spins in the limit m — > oo. In 
this case || we obtain a similar picture to that described 
above: namely a high temperature phase separated from 
a spin glass phase at low temperature. The form of the 
spin glass phase is RS with q non-zero and the equations 
for Qabcd may be solved to find that Q a bab stays critical 
below the transition with exponent 9 given by 8 = d — 1 
for RKKY-like interactions fUjfl , Applying the method 
to to = 1 Ising spins will lead to errors, but we hope 
that certain features will be correct. A case that can 
be analysed rigorously is that of a purely ferromagnetic 
interaction and Ising spins where we can demonstrate 
that the spin glass and ferromagnetic transitions must 
be simultaneous The Gaussian variational method 
fails in this respect, predicting a spin glass transition at 
slightly higher temperature than the ferromagnetic tran- 
sition. Indeed this effect has been noticed before, and 
Sherrington [T^ ] has identified relevant diagrams that are 
ignored in the GV approach. 

Another shortcoming, not related to the Gaussian vari- 
ational approximation, may also be present in our treat- 
ment of spin glass ordering. That is that we have only 
taken our analysis as far as order parameters with two 
replica indices which is known, for example in the Viana 
Bray model, not to be correct fTq] . This effect may be 
apparent in the case of an antifcrromagnetic interaction. 
There is no difficulty of principle in extending our meth- 
ods to consider operators with more replica indices, but 
in practice, the calculations soon become unwieldly. 

We would like to acknowledge useful discussions 
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FIG. 1. Phase diagram for the 3 dimensional RKKY in- 
teraction using Hartree approximation. 
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